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Abstract 

I In the framework of the A:T-factorization approach we have calculated in the frag- 

' mentation model the p^-spectra of Be mesons at the energies of the Tevatron and the 

. LHC Colliders and at the large pT domain. We compare the obtained results with the 

Q I existing experimental data and with the predictions obtained in the collinear parton 

^ • model. 

Q: 
o ; 

^ 1 Introduction 

^ I Recently, Be mesons were observed experimentally These particles have a special 

^ ' place in the doubly heavy meson family, because they consist of quarks of different flavors. 

^ . The Be meson spectroscopy has the same features as the spectroscopies of charmonia and 

^ ! bottomonia, however the decay and production properties of Be mesons are unique. There 

^ I is no annihilation channel for Be mesons into hadrons via the gluons and, consequently, the 

O I lifetime of the ground state of the be system is large, given by r ^ 0.46 ps [2]. The main 

^ ' channel of Be meson decay occurs through the weak decays of the 6-quark or the c-quark, 

which are respectively about 25% and 65% of the total decay width. 
Oh. The Be meson hadroproduction cross section is approximately 100 times smaller than the 

r-| ! (66)-pair production cross section, because the two pairs of the heavy quarks are produced 

^1 in a parton subprocess. The mass spectrum and the decay properties of the Be mesons 

are studied in detail in the framework of potential quark models |3j, in the QCD sum rules 
methods |3] and in the OPE approach jH]. The predictions at leading order in for the 
Be meson production rates in the 77, 7p and pp interactions were obtained in the collinear 
parton model [3 H 1101 • 

The Be meson production cross section became large enough to study experimental (~ 1 
nb) only in the region of very high energy or in the region of very small x ~ Mj ^ ~ 10~^, 
where x is the argument of the gluon distribution function in a proton (the quark contribution 
to the hadronic production of Be mesons is very small). In this region the /cr-factorization 
approach [TTl |T2| is more adequate for describing the perturbative evolution of the gluon 
distribution function, which satisfies the BFKL [12] or CCFM evolution equations, in 
comparison with the collinear parton model, which is based on the DGLAP equation jl5j . 

In this Latter we have calculated Be meson p-p-spectra at the energy range of the Teva- 
tron and the LHC Colliders in the fragmentation model and in the framework of the fcr- 
factorization approach. In the fragmentation model we have taken into consideration only 
the main contribution originating from fragmentation of a b-quark into a Be meson [T^ . 
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2 The factorization approach 



In the /cT^-factorization approach ^TJ which generahzes the colhnear parton model to 
the region of small x, the hadronic and partonic cross sections are related as follows: 

where $(a;, |/crP, /^^) is the unintegrated gluon distribution function in a proton (unintegrated 
refers to the transverse momentum), l^j^rP is the virtuality of the initial reggeized gluon, 
ki = XiPi + ki^T is the 4-momentum of the initial gluon, ki^T = (0, fcj,T, 0) is the transverse 
momentum of the initial gluon, ipi is the angle between the gluon transverse momentum and 
the fixed axis OX in the plane XOY, s = xiX2S—\kiT\'^—\k2T\'^. In the numerical calculations 
we have used the following parameterizations for the unintegrated gluon distribution function 
in a proton: JB by Blumlein JS by Jung and Salam and KMR by Kimber, Martin 
and Ryskin [T^ . 

To calculate the amplitude of partonic processes in the /cj^-factorization approach the 
polarization vector for the reggeized gluon in the initial state is taken to be fn\ IT^ 

e'ikr) = ^. (2) 

The squared amplitude for one of the processes in Be meson production in the k^- 
factorization approach, namely 

g* + g*^b + b (3) 

was obtained earlier in JT]. Here we use our original results from Ref. 



\Mig* + g*^Q + Q)\^ = Sn^aliMn + M22 + M33 + 2Mu + 2M13 + 2M23), (4) 

where 



Mn = ^('^i(4/3i - + 4- 4/52«;2) + 4A/?2(«;i - 2Pi){4kI + i) + 4/€?(A/t2 - Pi) 

+ AI31{kI + A(3iKi - 4A/€i/€2) + 4Ak2(A/€2 - 2/3i)(/s:2 - A^l) - 4/3^(2/52 - ^2)' + 4 
+ 4A/t2(2A/?2 + A)(2/t2 + ^) + 8A/3i/t^ + 4A/€2(S + u)(Ki + A/t2) (5) 



Moo = 



(4 - AP24 - ^I3l{4 + 4) + 16/3i/32'(/€i - (3i) + A\{Al3^4{n2 - 2P2) + iu 



+ Ki{4 



2P24 + {^^2 - /?2)(s + i-APl) - P2U) - A (2/32 - ^^2)u) - \^\^0i4 

+ 4A'/€i(4A/€2 - Kx{4 + 7i) + 2/3i7i) + K2(4/32 - «;2)((2/3i' - K^f + 4) + 44pi)i,^) 



2 



- A^(ki + 4 - Sfh^l - 8/32^2 -s'^+P + 8(3iKi{kI - 4(32^2 + i-u)- 2iu + 

- 2k\{kI - AP2K2 - 8pl + i-u)+ 2k^(8/?2 + t - u) - 8/52^2 - u)) 

- 2\{S(31kiK2 + 2(32{nl - A(3ikI - A(3ikI - + i-u)) + ^1^2(8^1 - s))) , (7) 



M12 = ^(iPiKl + Kl{s-8Pl-i-u) + 32PMPiK2-M2)+s'' -P -u" 

4:8tU V 

- kI{kI - n\- 8/32^2 + 8(5l - 2s) + 4/3i/ti(2/t2 - 8/32^2 + 8^1 - k\) + 2\^{k\ 

+ kI- s"^ + F + + AI32K2U + 2Kli + 2Kl{4:p2H2 +U)+ 4:j3iKi{2Kl - 4:p2K2 

+ i)j - 4x(ApiK2iKi - 2/3i) - K2(-4/3i«i + Pi{^kI + kI-s)+ kis) - /32(«? 

- ApiKl + Ki{3Kl + 8(3l-s)+2Pi{s-3Kl-Kl))jy (8) 



Mi3 = ^ f 2(A^ - l)^^ + 2k?(2/3i - Xk2 - S/^iA^ + 3X^2) + 2/3iA(2k^ + 4^2^ 

+ 4/32(t - - K2{16pi -2u))- kI{2{1 + A2)fi:2 - 2((5A2 - 2)p2 + GPiX)k2 

+ 24A/3i/32 - 8P^ + 2s- A^s + XH +2u- SX^u) + 2k2/92(8/3i +t + u)- 4(31^2 

- H2U + A^(4«;2 - 22/52^2 - + (£ - w)^ - K2(s - 32/^2 - 5f + 3w) + 2/32^2(5 

- lt + bu)) + 2ni{(3iX^{bnl - 16/32^2 - 3s + f - 3?2) + A (4^2 - 8/32 + 2s + f 

+ m) + A(3fi;^ - ll/32fi;2 + ^2(16/32^ - 8/3^^ + 2t-u)+ (52{lQ(5l + s - 3£+ 3m)))) ,(9) 

M23 = -77^ fA^(4Kf + 24 - 22(3ik\ - s^ + 2/3iKi(5«;^ - 16/32^2 + s + 5f - Til) 

- Ki(2k2 - 16/32^2 - 32/3i + s + 3f- 5u) + kI{s + 3i - tt) + - -u)^ + 8/32^2^1) 

- 2{kI{kI - 3/32^2 + 4/32^ + i) + K2(«2 - (^'^'4 + '«2(s - 4/3^2 + i) + /32(8/3i2 - s)) 

- /3iKi(8/32' - 2kI + i + £t)) + 2A(k?(k2 + 2/32) + /3ik?(4/32 - IIK2) + /3i(4/32(ii - £) 
+ 3^2 - I2/32K2 + «2(16/32 + s + 3£ - 3ii)) - Ki(3«;2 - 8/32^2 - «2(16/3i - 80^ - 2s 

- 3t + «)+/32(16/32-t)))), (10) 

where 



t = 




(11) 


■u — 


u-rriQ, 


(12) 


A = 


\Pt\ COS</?i, 


(13) 


/32 = 


|Pt|cos (p2, 


(14) 


Ki = 


\kiT\, 


(15) 


K2 = 




(16) 


X = 


cos {(fi - (f2) 


(17) 
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and pt — transverse momentum of anti-quark. Since these are more suitable for numerical 
calculations. In those calculations, which were in the framework of the collinear parton 
model, we used the GRV |21j parameterization for the gluon distribution function. 



3 The fragmentation model 

The analysis of the B,, meson gluon-gluon production in the framework of the collinear parton 
model shows [El 13 Ej the dominant role of the fusion mechanism up to px ~ 30 GeV. It takes 
into account the total gauge invariant set of the diagrams which describes the parton process 

g + g ^ B^ + h + c. (18) 

Thus, the fragmentation approximation, incorporating factorization of the heavy quark pro- 
duction process and the soft process of a meson created in the final state, applies only at 
Pt > 30 GeV. Of course, at the Tevatron Collider this region of pt is far from being available 
for experimental study. On the other hand, at the LHC Collider the pt < 50 GeV region 
will be under consideration. 

In the fragmentation model, the B^ meson production cross section can be presented as 
follows: 

da{pp BcX) = / dzDi_,Bc{z) ■ da{pp i), (19) 

where the sum is taken over all the relevant partons i = c, b and g (but the main contribution 
comes from 6-quark fragmentation). 

The fragmentation functions for the 6-quark splitting into a vector B* and pseudoscalar 
Be meson at the starting scale of the perturbative QCD evolution were obtained in Ref. 

[221 ES]: 

D,^sMh) = '^"''"'f (6 - 18(1 - 2r). + (21 - 7,r + 68r^).^ 

81m^ (1 — (1 — r)zj° \ 

- 2(l-r)(6-19r + 18r2)2^ + 3(l-r)2(l-2r + 2r2)/), (20) 

- 2(1 -r)(4-r + 2r=^)^3 + (1 -r)2(3-2r + 2r2)z4j, (21) 

where r = ^ ' l^-BclO)!^ is the squared Be meson wave function at the origin, which 
can be calculated using the nonrelativistic potential quark model |21]. In the numerical 



calculations, we used the result \E'b^(0) = y-^fscj where fs^ = 490 — 560 MeV is the Be 
meson leptonic decay constant. 

The QCD evolution of the fragmentation functions Di^^^ and D^^^^ are described by 
the DGLAP evolution equation 

dD . r,. a;o(u^) dx ^ ,x 

— {Z,fi') = ^^ / — 

a/x^ Zn J X 

where Pq^q{z) is the standard quark-quark splitting function. 
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4 The results 



The results of our calculation for the p-p-spectra of the B* meson are shown in Fig. ^ for 
the energy of the Tevatron Collider ^/s = 1.8 TeV, and in the Fig. El for the energy of the 
LHC Collider y/s = 14 TeV. The spectra that we obtain are compared with those obtained 
earlier in the coUinear parton model via the fragmentation mechanism. The curves in Fig. 
n] correspond to a choice of the parameters from Ref. = 1.5 GeV, = 4.9 GeV, 

as = dsiPT + "^bJ Ibc = 490 MeV. The curves in Fig. |21were obtained using the set of 
parameters from Ref. 0: = 1.5 GeV, rrib = 5.1 MeV, as — 0.23 and fs^ = 560 MeV. 

Fig. n] shows that our result obtained in the coUinear parton model agrees well with the 
result obtained earlier in Ref. pH] • The curves obtained in the fcr-factorization approach lie 
approximately 2 — 5 times higher than the collinear parton model prediction in the region 
Pt > 10 GeV. The maximum value of the Be meson production cross section is obtained 
using the JS [TH] parameterization of the unintegrated gluon distribution function, and the 
minimum value is obtained using the JB J7] parameterization. 

Note that the slopes of the p^-spectra obtained in the collinear parton model and in the 
/c-r-factorization approach are equal, and the situation is the same as in the case of D* meson 
production in ep interactions ^251 OB] . 

At present, there is no experimental information on the pT-spectra of Be mesons. It is 
known that the integrated production cross section in the region pt > Q GeV, |?/| < 1 and 
^/s = 1.8 TeV is given by 0"^^ ~ 10 ± 6 nb As shown in Ref. [H], the uncertainties in 
the calculations are about ~ 50%. This value depends on the choice of the masses, the 
constant and the leptonic decay constant Jb^- Furthermore, our calculations performed in 
the /cT^-factorization approach show that the B^ meson production cross section has a strong 
dependence on the choice of the unintegrated gluon distribution function — a variation by 
a factor of 2 was observed. 

The obtained values of the integrated Be and B* meson production cross section at the 
energies of the Tevatron and the LHC Colliders are presented in Table ^ One can see 
that, with the fragmentation model, the results obtained in the collinear parton model are 
smaller to those obtained in the fcy-factorization approach. The theoretical predictions are 
shown together with their uncertainties, which are significantly larger in the /cr-factorization 
approach. The absolute value of the cross section obtained in the collinear parton model 
{asc = 1-7 ± 0.8 nb) is much smaller than the CDF experimental data {aB^ = 10 ± 6 nb) 
PQ 12] . On the other hand, the prediction of the /cT-factorization approach agrees well with 
the data {cb^ = 7.4 ± 5.4 nb). 

Taking into account the relative roles of the fusion and fragmentation mechanisms in 
Be meson hadroproduction, one can suppose that the cross section calculated in the fusion 
model and in the /c-r-factorization approach will be more large than that obtained in the 
fragmentation model and in the /cT-factorization approach. The calculation of the Be meson 
hadroproduction rates via the partonic subprocess p8|) with the reggeized initial gluons will 
be presented in our forthcoming paper. 
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Table 1: The summed Be and B* meson production cross sections in the different models at 
the energies of the Tevatron and LHC Colliders. The data from the CDF Collaboration PflEj 
are shown. The cross sections are in nb. 





\y\ 


PT,miny 


Parton 


/cT-factorization 


Experimental 


TeV 




GeV 


model 


approach 


data 


1.8 


< 1 


6 


1.7 ±0.8 


7.4 ± 5.4 


10 ±6 


14 


< 2.5 


10 


28 ± 14 


122 ± 90 
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Figure 1: The B* meson pr-spectrum at ^/s = 1.8 TeV and \y\ < 1 in the fragmentation 
model. The continuous hues are a result of calculations in the /cT-factorization approach 
with the different parameterizations of the unintegrated gluon distribution function. The 
dashed line is the prediction from the coUinear parton model. The points show the results 
obtained in the Ref. [TEj . 
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Figure 2: The B* meson pT-spectrum at ^/s — 14 TeV and \y\ < 2.5 in the fragmentation 
model. The continuous hues are results of calculations in the fc^^-factorization approach with 
the different parameterizations of the unintegrated gluon distribution function. The dashed 
line is the prediction from the coUinear parton model. 
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